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INTRODUCT | ON 


The princlpal objective of these programs when first written was to 
reduce marine geophysical data from thelr raw form (readings of a strip chart 
recorder or Precision Depth Recorder) and combine them with ships’ naviga= 
tlon Information to obtaln fully reduced values of the geophysical para= 
meters (free-air gravity, residual magnetic Intensity or depth) and the 
position at which the values were obtained, Subsequently the system grew to 
Incorporate methods for computer adjustment of ships! track, for plotting 
the values as profiles of various kinds, for annotating the values on Merca= 
tor charts, for averaging values over various areas, etc. An extremely 
Important function of our present system Is the storage of large amounts of 
geophysical data and the abl itty to retrieve them quickly and efficiently 
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from the random access disk cartridges, 

Our present computer Installation consists of an IBM 1130 (with a 
basic cycle time of 2,2 microseconds and a storage capacity of I6k 
l6=bIt words)which tncludes 

(1) Two disk drives which accept Interchangeable disk cartridges 

(11) Printer, with speed up to 350 IInes per minute 

(1'1) Card reader, with speed up to 1000 cards per minute 
(1V) Card punch, with speed up to 120 cards per minute and 
(V) 30" wide Calcomp digital Incremental plotter which plots at 
200 steps/second, each step belng !/100 Inch, 

All the programs described here are In Fortran 2 and have been comp! led 
and run with the IBM 1130 Monitor 2 system, 

Thus we have a medium size and speed computer Installation with on-line 
Input and output from cards and disk, Our cholce of this Installation was 
gulded by two prime considerations, Firstly, for a majority of the Jobs the 
output Is a graphical plot. To have a very fast computer would mean that 
for a large part of the time the computer would be walting for the plotter, 
a procedure that Is wasteful, It ts, of course, feasible to use a slower 
computer than the 2,2 microsecond cycle time IBM 1130, In fact, we had 
previously used the IBM 1620 MK |, the IBM 1620 MK I! and then the 3,6 micro= 
second cycle time IBM 1130, The change to a faster IBM 1130 was made to 
speed up data reduction to an optimum level for graphical and analytic work, 

Secondly, we find that the Interchangeable disk cartridges with random 
access are quite sulted to our needs, Magnetic tape Is, of course, cheaper 
In terms of the amount of Information stored, but random access Is not 
possible with It. NoneInterchangeable disks or drums allow for faster access 


times but thelr total storage capacity Is |IImited and they tend to be more 
® 


expensive than Interchangeable disks, 


An alternative scheme to ours Is one that would ut!!!ze a computer 
with a much higher speed and with output on magnetic tape. An off.l|ine 
plotter would be used for plotting, We have not utIl|!zed such a system 
for three reasons, One, the cost of an off=|Ine plotter system Is not 
appreciably less than that of our computer Installation, and It Its readily 
seen that a computer with an on-l!ne plotter Is much more useful than an 
offline plotter with only a plotter control, Secondly, as wlll! be seen 
later, we use the computer In many cases as a ||nk between a compact form 
of data storage and display on a plotter, A very high speed computer which 
would read data In a certaln format from magnetic tape and then essentially 
write It In a different format on magnetic tape for controlling the plotter 
will not tncrease throughput. Thirdly, In our case, output on the plotter 
Is used also In editing and error detection during data processing. Rapid 
turn-around times are achleved and wastage of machine time Is ellminated by 
having a medium speed computer with an onwl!ne plotter devoted almost 
entirely to reduction of marine data, In our opinion, a very high speed 
computer wlll become useful for processing data such as ours only when a very 
effictent method of time sharing Involving large amounts of data transfer as 
well as an efficient (for a fast computer) plotter for making plots of large 
dimensions (30" x 100') on=lIne, Is developed, 

The |nput data for our programs have generally been hand digitized, In 
some cases we have digit!ized analog records with an automatic XY digitizer, 
There Is no reason why raw data obtalned tn digital form by a shipboard data= 
logger could not be ut!lIzed, 


Agatn, while we have used the present system In the laboratory, there Is 


no reason why It could not be used on board a ship. In fact we have used a 
few of the programs described here for a smaller computer (PDP=8S) and 

used these on board R/V VEMA and R/V ROBERT D. CONRAD (Talwan!, Dorman, 

and Kittredge, 1969), WIth some modifications, our system could be made 
Into a real time data acquisition, reduction and display system on board 
ship. However, we do not favor real time use of computers on board ships 
for data reduction because, among other factors, we feel that 


(a) Unless adequate tlme=sharing procedures are developed, the 
computer will not be utllIzed efficiently 


(b) The necessity for editing of data precludes real-time processIng 


(c) Adjusted navigation essential for the final data reduction can not 
at present be attained In real time 


(d) The system described here can be used on board shIp as soon as 
navigation data are avallable to complete adjustment of track, 
which means that final reduced data are available after a few hours 
of data acquisition, If destred, 

It 1s felt that the most Important way to speed up reduction of underway 
martne data Is to develop good digital data logging systems, The data reduc» 
tion can be done elther a few hours later on board ship or a few weeks later 
In the laboratory, These remarks apply to data with very low sampling 
rates = a data polnt every minute or less often, Where data rates of several 


hundred numbers per second are Involved, on |Ine computers with real-time 


processing may well represent the most efficltent method of data handling, 


NAVIGATION POINTS 
Basic to a large number of programs described here are the navigation 
data represented by navigation points. A navigation potnt Is defined as a 
point where a change In speed and/or course of the ship occurs; In other words, 


the ship's track can be represented by segments at constant course and speed 


between navigation points. Stnce the ship's speed or course Is never 
absolutely constant, one must state how large a change of speed or course 
Is needed to define a navigation polnt. We have adopted the convention 
that whenever the ship actively maneuvers, that !s, the engine speed or 
the course steered Is varied, a navigation polnt results, In addition, 
track adjustment at navigation fixes may result In apparent course and/or 
speed changes, These are also considered navigation points, Examination 
of distance covered by the ship through water (by an e.m, log) Indicates 
changes of speed, When the change of speed (through water) exceeds about 
O0,! knot, a navigation polnt Its assumed, Finally, our method of Inter= 
polation between navigation polnts makes a flat earth approximation and the 
navigation points must not be very far apart, Hence additional navigation 
points are picked so that the distance between them does not exceed about 


30 nautical miles, 


SIMULTANEOUS INPUT OF NAVIGATION AND GEOPHYSICAL DATA 
A large number of our programs uti|l!ze navigation data together with 


raw digitized geophysical data (gravity, topography, or magnetic data) 


presented as a time series, In all our programs the data are allowed to be 
at unequal time Intervals and there Is no requirement that the different data 


all be read at the same times, The procedure that one would Intuitively 
follow would be to have these Input data on cards (paper tape or magnetic 
tape); to load these Into memory; and then to perform the data reduction. 
Such a procedure calls for a large memory storage and while thls Is no prob= 
lem for the very large computers avallable now, the smal! computers that we 
used Inttlally did not have large memory storage capacities, Since In the 


data reduction the geophysical data corresponding to a particular time 
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uti ltzed only the navigation data close to that time, It was possible to 
perform the data reduction without storing all the data In memory, Instead, 
two separate Input devices were used, one for navigation data, the other 
for the geophysical data, A few geophysical data polnts were read Into the 
computer core memory together with the navigation polnts required; the two 
Inputs were used In parallel and thus considerable savings In memory storage 
were reallzed, We have carrlted over thls same procedure to our present 
system where the two Inputs are stored on disk and then called for as neces= 
sary, It Is clear that many modtflcations of thls are possible for other 
computer systems, It Is polnted out, however, that If only one Input device 
(such as a stngle card reader) Is used, the navigation data will have to be 
stored first In core memory and hence there should be provision for an 
adequate memory storage capacity, 

The specific procedures for loading gravity, topographic, and magnetic 


data on disk are described In a later section, 


STORAGE OF DATA 

Ortginally our output data from the primary reduction programs cons!sted 
of one punched card per data polnt., On this card would be punched the latitude, 
longttude of the polnt where the value was obtalned together with the value, 
and such data as distance along track In nautical miles, As our output data 
accumulated, storage of cards became a problem and we began to look for ways 
to make the storage format more compact, It Is Important to reallze that the 
exact form of storage of output data Is not Important, What Is Important Is 
that the storage be compact and that It be easlly retrievable for programs 


which display It or use It In other ways, We therefore declded not to store 


the output data separately at all but to use the stored Input data In Its raw 
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form for use In all programs that would have ut!|!zed the output data, 


Compactness Is real!zed because the Input data are In a packed form, The 
machine time for the actual reduction of data is trivial so that no appre= 
clable time ts lost In reducing the data each time they are used for any 
purpose, Our stored data consist of two time serles. One serles consists 

of navigation points and gives the ship's posItton at various times, The 
other gives the value of the geophysical parameter In Its raw form at varlous 
times, together with the constants necessary to convert It to the reduced 
values of the geophysical parameter, Of course, the time serles giving the. 
geophysical parameter need not be In thelr raw form, They could be In thelr 
reduced form, The gravity programs, for Instance, can use Input data In the 
form of free=alr anomalles rather than chart readings, The magnetic programs 
can use {nput data In gammas, The reason why we do not usually keep the data 
In a reduced compact form Is that then we would have to store two sets of 
data (raw and reduced), whereas now we require storage for Just one set of 
data = the Input data for use In all our programs, 

Some of our programs can also use output data punched one polnt to a card, 
These data are generally referred to as answer cards or output cards on the 
comment cards which describe the programs, 

Some specific detalls of the organization of our data storage are pere 
haps of Interest. The storage of data !s by ships’ "legs", A three to five 
week trip from one port to the next constitutes a leg. Two data filles each 
are used for the storage of navigation, gravity, topography, and magnetic data, 
For VEMA Crulse 24 leg 13, these flles are named A2413, 82415 for navigation 
data; C2413, 02413 for gravity data; E2413, F2413 for topographic data; 

G2413, H2413 for magnetic data, The exact contents of the files are described 


later, Take, for Instance, an area such as the South Atlantic, Three disk 
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cartridges (let us call them #2006, #3006, #4006) contaln all Lamont data 

In the South Atlantic. All the navigation filles for all the legs of 

various crulses are loaded on disk cartridge #2006, All the gravity and 
topography data filles are stored on disk #3006, All the magnetics data 

flles are loaded on disk #4006, All the data reduction programs are loaded 
on disk #2006, Our computer system has two disk drives, By mounting carte 
ridge #2006 on drive 41 (logical drive #0) and #3006 on drive 2 (logical drive 
#1) all the gravity and topography data In the South Atlantic are avallable 


as Input Into all the programs described here, 


NAVIGATION PROGRAMS AND LOADING NAVIGATION DATA ON DISK 

As mentioned earller, navigation data consisting of navigation polnts 
are used as Input to a large number of data reduction programs, The naviga- 
tlon poltnts are stored on two filles (for definition of Flle see IBM 1130 
manual; also see Table |) on disk elther by the DSKNV or the FIXPI programs, 
The data stored are basically the time and pos!Ition at each navigation point, 
Time {s stored tn hours from January 1 of a base year which can arbItrarlly 
be defined for any leg. (Base year Is defined by a quantity JBSYR that Is 
read In, Time Its reckoned from 00 hrs G.M.T. of January 1 of the year JBSYRt1, 
If no base year data are read In, JBSYR Its set#1959,) Latitude and longitude 
are stored In degrees and fractions, The time zone Is also stored so that 
ship's tlme can be recovered, If FIXPI program has been used, the speed and 
direction of ship's drift obtalned In the track adjustment procedure are com= 
bined Into one word (azimuth of drift In degrees times ten plus speed of 
drift) and also stored, DSKNV program Is used when the navigation polnts are 
avallable from an adjusted track, The Input to DSKNV Is the time and position 
of navigation potnts, FIXPI Its used when track adjustment Is performed by the 
computer, The Input data In this case are "Pitlog polnts" (which are the 
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times and pitlog, or e.m, log, distances of polnts where the ship changed 
course and/or speed together with the course steered subsequent to the 
time of the "Pitlog Point") and flx data (which consist of the time and 
position of fixes obtalned by the shlip by any method whatsoever), FIXPI 
uses the usual method of track adjustment = a dead reckoned track Is ob= 
talned from one fix to the tlme of the next flix, The difference between 
the position at the second flx and the dead reckoned position at thls time 
Is attributed to a constant "current" between the times of the two fixes, 
The constant "current" {s applied to the navigation potnts between the two 
fixes to obtaln the adjusted track, 

The MGNAV program (Iisted !n Magnetics Programs) Is used to determine 
the value of the reglonal fleld as given by Caln et al, (1964) at each navie 
gation point and to store {tt on the disk fille, The purpose of doing so Is 
to be able to obtaln the regional fleld at each magnetics data polnt by Inter- 
polation between navigation points, rather than by an Individual calculation 
at each data polnt. MGNAV programs require the polynomial coeffictents to be 
avallable on disk as a stored file. This fille CF1 must therefore be stored 
previously, which Its done by a program CCF1, Instead of using the reference 
fleld of Catn et al, (1964), It Is possible to use the new International Geo= 
magnetic Reference Fleld adopted In 1968 by usIng program NIGRF In place of 
MGNAV and ustng a different set of polynomial coefficients to load as file 
CF2, 

NVANS program ts a check for proper loading of navigation data on disk, 
It takes Its Input from the navigation file on disk, The output gives the 
time, latitude, and longitude of navigation polnts, distance In nautical miles 
along the track, speeds and courses between them, and If desired, the regional 


magnetic fleld at each navigation polnt. 
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ETVOS program also uses the stored navtgation data on disk and outputs 
the Eotvos correction (used In gravity reduction) on the preliminary plot 
(a plot of gravity and topography used during data reduction for checkIng 
purposes), 

Programs PLODI and PLTIM also use the navigation data stored on disk 
and are used to make the distance (at the bottom) and time (at the top) 
annotations on the final plot (the plot on which reduced values of gravity, 
topography, and magnetics are displayed In proflle form = see Fig. !). 

Program MERN2 (with I!nk to MERN3) Is used to plot the ship's track 
on a mercator chart at any destred scale. Varlous options with regard to 
annotations of hour marks, distance, and time are avallable, The use of 
MERN2 ts demonstrated by two examples, The bathymetric contours In the 
second example were obtalned from DIGM3, a program that Is not described 


In this report, 


LOADING GRAVITY, TOPOGRAPHIC, AND MAGNETIC DATA ON DISK 

For each leg of the shlp we store separately the gravity, topography, 
and magnetic data, The topography data, for example, consIst of a time 
serles that defines POR readings at various times (not necessarl ly equal ly 
spaced), To take advantage of the fact that some quantities In the time 
serles do not change rapidly, we have Introduced "day cards", Punched on 
the topography day cards are: Day, Month, Year, Time Zone, number of topo= 
graphy data cards (LB) that day, number of polnts (LC) on the last card, 
Matthews’ Zone number (NMZ) for the ship that day. Followltng the day card 
are a number of topography data cards for that day, which In a packed format 
have time (hr., min.) and depth (nominal fathoms =1/400 sec, reflection 


time) punched on them, This Is followed by the second day card with Its 
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data cards and so on, A blank card terminates the deck, Program MTOPD 
loads this deck of topography Input data; the day cards on to an "E" flle 
(one record per day card = Index KS) and the data cards on to an "F" #lle 
(one record per data card = Index KT), By "E" flle we mean that for leg 
#13 of Crulse 24, the convention adopted Is to name the flle E2413, SiImle 
larly, the "F" fille Is named F2413, Stnce the two files are kept separately 
In storage, It Is necessary to store some Information on the day cards file 
to determine which records on the data cards fille are on that day. If KT a 
Is the record number of the first data card on a certain day (m) the 

number KT,=1 ts also stored as a part of that day record KS,,, To monitor 
that the data have been stored properly MTOPD prints out a |Isting of the 
day record file, that Is Day, Month, Year, TZ, LB, LC, NMZ, KT,=1, KS for 
each day, 

MATZ Is used to store the velocity of sound In the various Matthews! 
zones on disk. This flle has, of course, to be stored only once on each disk = 
and not stored separately for each leg, 

MGRAD similarly stores gravity data onto disk on "C" files for day cards and 
"D" f#lles for the data points, The data potnt time serles consists of time 
(hr., min.) and strip chart readings, Constants JDIAL, DRR, CF, GZ are 
needed to convert the strip chart reading JG(I) to a gravity value GRAV by 
the formula: 

GRAV=GZ + (JDIAL + JG(I) x DRR) x CF x 0,001 


The constants are used for Graf=Askanla meter data reductions, and are 


stored as part of the "gravity day" card data, If any of the constants change, 
a new "gravity day" Is defined, Thus, In general, there are more "gravity 
days" than actual calendar days, 


As mentioned earller, we store our Input data not only for data reduction, 


but also for final data display, plotting, etc. If one starts with reduced 
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data (gravity values or freesalr anomalles) rather than strip chart 
readings, the program MGRAD Is used equally well to store the reduced 
data. In this case, the constants JDIAL, DRR, CF, GZ are omitted or have 
different meanings, The programs that employ the stored data have pro= 
vistons to use constants appropriately, 

MMAGD, similarly, stores magnetic data, Again, the "magnetic day" 
cards store the constants used In the data reduction and a new day card Is 
used whenever any of these constants change, MMAGD Is used regardless of 
whether the Input data ts In gammas or whether the time T for a number of JCY 
cycles of the proton precession Is measured, In the latter case, the fre= 
quency !s JCY/T, The time T Its measured by counting the number of cycles N for 
the same time Interval of a standard frequency source (let frequency be JOSC 
In untts of hundreds of KC), that Is, T = N/(JOSC x 100,000), The two high 
order digits of the count N are called JCT and stored with the "magnetic 
day" data; the three low order digits C(J) are stored as the data time 
series, Thus, total Intensity fs 

H = (JOSC x 100,000 x JCY x 23,48738)/(C(J) + JCT x 1000) 

JOSC and JCY also are stored In "magnetic day card records," 

If the Input data are In gammas, the two high order digits are termed 
JCT and stored In the "day records" and the three low order digits become 
C(J) and stored as the data time serles, Provision has to be made In the 
programs that use the stored magnetic data to discriminate between data 
stored as perlod or as gammas, 


Some detalls regarding the data storage flles are given tn Table |, 
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BASIC GRAVITY, TOPOGRAPHY, MAGNETIC DATA REDUCTION PROGRAMS 

The Input for these programs consists of the navigation data and 
the gravity or topography or magnetic data stored on disk, The output 
consists of printed reduced values and plots at varfous scales, The 
printouts are for gravity (GRDSK) day, month, hours and minutes, latitude, 
longitude, distance along track In miles, distance along track In km, 
observed gravity, International formula, free=alr anomaly, crulse Iden-= 
tification, and Eotvos correction, For topography (TPDSK) depth In nominal 
fathoms, corrected fathoms, and corrected meters and the year are printed 
Instead of gravity data, For magnetics PRMG2 prints out hour and minute, day, 
month, time zone, latitude, longitude, distance along track In nautical 
miles, observed magnetic total Intensity, theoretical magnetic regional 
value of field, residual anomaly, last two digits of the year, and crulse 
Identification, 

As a part of our standard data reduction procedures we have found It 
useful to plot the reduced data on preliminary plot form at the fixed scales 
of 50 nautical miles to the Inch and 50 milllgals; 500 fathoms; or 500 gammas 
to the Inch, This Is also accomplished by the programs GROSK and TPDSK, 
and PRMG2, Long profiles are automatically broken up fnto two parts (this 
has been called folding) and plotted on top of each other, The zero for the 
free-air gravity plot lies, for Instance, at 8 Inches from the bottom of the 
plot and at 23 Inches for the "folded" plot, As mentioned earller, ETVOS 
plots the Eotvos correction on the preliminary plot scale, The plots on the 
same paper of free-air gravity, topography and Eotvos are extremely useful In 
checking both the gravity and topography data for errors, We generally find 
that two or three preliminary plots are essential before the data are con- 
sidered reduced In final form, 

Final plots are plotted on specially printed paper by programs GRDSK, 


TPDSK, and MGFN2, Such a plot Is shown In Fig. |. Navigation data (miles 
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Figure | 


Final plot of Depth (flnest I!ne), Freesalr Gravity and Residual Magnetic 
Intensity (thickest IIne) obtained automatically by the on=IIne plotter, 
ustng programs PLODI, PLTIM, GRDSK, TPDSK, and MGFN2, Orlginal scale of 
the plot fs 4" = 20 milligal = 200 nominal fathoms = {00 gammas = 20 
nautical miles, Scales for Gravity, Depth, and Magnetics are marked G, D, 
and M Respectively, Program PLODI and PLTIM describe the other annotations, 
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and time, etc.) are also plotted on the final plot by programs PLODI 

and PLTIM, The scales are always kept the same for the final plots; 
however, the zeros for gravity and magnetics are shifted to avold con= 
fuslon between Individual traces, (Prefinal plots are plotted by programs 
GRDSK, TPDSK, and PRMG2; these plots are merely scaled down versions of 
the final plot and are used to determine the optimum zero shifts for 
gravity and magnetics.) 

The residual magnetic anomaly Is used In all these plots, It Is 
computed by Interpolating a regional value between values computed for the 
regional fleld at navigation polnts and stored on the navigation file by 
the program MGNAV (or NIGRF), Two polnts should be noted, It Is obvious 
that the regional magnetic fleld flle must be estab||!shed by MGNAV (or NIGRF) 
before the PRMG2 or MGFN2 can be used, Secondly, the procedure described 
above achieves considerable savings of time at trivial sacrifice of accuracy 
Ca gamma) over other methods In which the regional value of the magnetic 
fleld ts computed at each data polnt, 

BOUG uses the Input gravity and topography data disk flles and computes 
the simple Bouguer anomaly at all polnts where gravity values have been 
obtatned, The printed output Includes the free=alr and Bouguer anomalles and 
the corrected depth In meters as well as the depth In nominal fathoms at the 
polnts where gravity data were obtained, 

Note that for all the data reduction programs the positions for the 


geophysical data points are obtained by Interpolation between navigation 


points; hence there must be no geophysical data points earlier than the first 
navigation polnt or later than the last navigation polnt. 


Subroutine DMAG2 Is used In all the magnetic data reduction programs, 
Its function ts to read the navigation and the magnetic data and do the 
magnetic reduction to obtain the position and the residual and total Intensity 
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values, The corresponding functions In the gravity and topography 
programs have not been separated out Into subroutines (as they wel! might 


have been) but are repeated explicitly In every program, 


PROGRAMS FOR ANNOTATION OF GRAVITY, TOPOGRAPHY, OR MAGNETIC VALUES ON 
MERCATOR CHARTS 


Programs MGRA4 (with IInk to MGRA5), METC8 (with I!Ink to METC9), and 
MGGG8 (with I!nk to MGGG9) respectively, annotate gravity, topography, or 
magnetic values on Mercator charts, The Inputs are Identical to the 
Inputs used for the basic data reduction programs, The outputs are anno= 
tated values on the Mercator charts at any deslred scale, There are options 
In METC8 to annotate nominal fathoms or corrected meters and In MGGG8 to 
annotate residual (In gammas) or total Intensity magnetic anomalles (In tens 
of gammas), 

METC2 (with |Ink to METC3) Is a modification of the topography annota= 
tlon program METC8 for use In contourtng. A contour Interval Is specified; 
the position of every multiple of the contour Interval Its Interpolated and 
represented by a tick mark along the track, If the contour Interval Is 
50 fms, then the tick marks corresponding to every 500 fms are annotated, 
Annotations are also made where reversal of slope makes the tick marks ambl= 
guous, Again If, for example, the Interval Is 50 fms and the Interpolated 
depths are 1450, 1500, 1550, 1600, 1650, 1700, 1750, 1800, 1850, 1900, 1950, 
2000, and 2050 fms the regular annotations are at 1500 and 2000 fms and In 
addition the second time any Interpolated sounding Is repeated, that Is, the 
second 1650 Is also annotated, 

All the above programs use I|Inks, In each case the first of the | Inked 
programs merely draws the Mercator grids and stores the varlables necessary 


for the Mercator plots. The IInks are devices for saving storage, Note that 


21 


both of the I!nked programs call for Input cards, 
The subroutine MERGR Is used by all the annotation programs except 


METC2, 


PROGRAMS FOR PLOTTING OF GRAVITY, TOPOGRAPHY, OR MAGNETIC VALUES AS PROFILES 
ALONG SHIP'S TRACK ON MERCATOR CHARTS 


Programs MGRA6 (with I|Ink to MGRA7), METC6 (with Itnk to METC7) and 
MGGG2 (with IInk to MGGG3) respectively, plot profiles of freeealr gravity, 
depth and residual magnetic fleld Intensity along the ship's track, In 
plotting the profiles the following cholce Its avallable, Elther one can plot 
values along a direction perpendicular to the ship's track, or along a fixed 
azimuth, For example, let the ship's track run through 10N 5E tn a NE direc= 
tlon and the gravity value at this polnt be 100 mgal to be plotted at 100 
mgals to the Inch, If the first alternative Is taken, a polnt corresponding 
to the gravity value will be plotted 1" NW of 10N 5E on the Mercator chart, 
lf the second choice Is taken, the point Is plotted at 1" along azimuth THET 
(read In as data) from {10N 5E, 

The Inputs for these programs are Identical to the Inputs used for the 
bastc data reduction programs, There Is an option In the METC6 program to 
use elther uncorrected fathoms or corrected meters for the plot, 

The subroutine MERGP Is used In the above programs, It has options 
which allow the plotting of the track Itself and Indication by crosses on the 


plot of the times at which navigation polnts occur, 


PROGRAMS FOR OBTAINING PROJECTED PROFILES OF GRAVITY, TOPOGRAPHY, 
OR MAGNETIC VALUES 


Programs NAGS1, NATS2, and MGPR2 respectively project values of free=air 
gravity, depth, and residual magnetic fleld Intensity along a |!ne which 


passes through a specified point at a specified azimuth, 
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Figure 2 


Flg, 2 (top) shows a Mercator chart with ship's track ABCB"D, The 
| tne on which one wants to project, say, the gravity values Is drawn through 
O (VLAT, VLON) with azimuth ANGLE. The gravity value g, at A Is projected 
to A' and so on, Note that there would be an overlap of projected values 
for ship's track BCB"; this Is avolded through programming and results In the 
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discrepancy between 9, and Ion In Fig, 2 (bottom), If one wanted to obtaln 
the profile corresponding to CD, for Instance, one would have to start the 
projection at C Instead of A by appropriately choos!ng XMIN (minimum dis= 
tance along track tn nautical miles at which projection starts), 

Where the ship's track and polnt O are within 1° of latitude of each 
other the flat earth approximation ts used, For other cases a correction 
Is applted (see TROJ2 subroutine), 

Projected values and projected distances (In km) are printed out and 
punched out on cards In addition to belng plotted, 

Subroutine TROJ2 Is used fn programs NAGS1, NATS2 and MGPR2, The Inputs 
for these programs are Identical to the Inputs used for the basic data reduc= 


tlon programs, 


PROGRAM FOR OBTAINING AND PLOTTING AVERAGE FREE=AIR GRAVITY VALUES 

AVFRE ts used to average freeealr gravity over given sized "squares" 
of latitude and longitude, Raw gravity data for a ship's leg stored on disk 
can be used as Input data to obtaln average values for that leg, The averages 
are printed out as well as punched out, The punched output cards can them= 
selves be used as Input to obtaln averages for different legs and for larger 
squares, 

AVPLT ts used for plotting the obtalned averaged free-air values (out= 


put of AVFRE) on a Mercator chart, 


SYSTEM SUBROUTINES AND FUNCTIONS 
Subroutines and functlons used by two or more programs are described 
here and are |!sted as system subroutines, Other subroutines are |Isted 


together with the corresponding maln programs, 


Subroutines for calculation of time, 

AMAY calls the tables of AM and AY from storage, These are 
respectively the number of accumulated days In the months of a year and the 
number of accumulated days In years starting from January 1, 1960, AMAY 
must be called prior to calling TXX or DXX. TXX computes time to 0000 
hours of given day, month, year from January 1 of JBSYR+1, For example, 
If JBSYR = 1963, times are calculated from January 1 of 1964, DXX operates 
In reverse fashion to TXX, If the time In hours (and fraction) from a 
known base year JBSYR+1 Is given, DXX calculates the day, month, year, 
hours, and minutes, Tlme zone Is needed to obtain ship's time rather than 
G.M.T. 

These subroutines are written for years 1960-1970, For years earlier 
than 1960, programs using TXX will still work; however, application of DXX 
to times obtalned by TXX for these years will yleld erroneous dates, For 


many purposes, this error Is Inconsequenttal, 


Subroutines for calculation of distance, course, and conversion of degrees 
and minutes Into fractlons and vice versa, 


COMDI computes the distance between polnt 1 (Latsy,, Long=x, ) and 


point 2 (Latey., Long=x.) by the formula 


((4e 4.) + (He %) cos CC 4% #) /2))" 


Note that this ts-a formula uti lIzing a flat earth approximation, 
Given a nerth vector XXN and an east vector XXE, CATAN computes In 
degrees thé azimuth of the resultant vector measured In clockwise sense from 


the north, 
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Given a coordinate tn degrees and minutes, DEGFR converts It Into 
degrees and fractions, Conversely, DEGMI converts a coordinate In 
degrees and fractions to degrees and minutes and fractions, 

DELON determines a difference In two longitudes so that the absolute 
value of the difference Is always less than 180°, to avold problems on 
crossing 180° meridian, 

Subroutines for Mercator and other plotting, 

MEGRI draws the outl!ne of a Mercator grid, glven WLONG, ELONG, SLAT, 
and ANLAT, the western, eastern, southern, and northern boundarles respec= 
tively at any specifled scale of SCALE Inches per degree longitude for any 
latitude range between 81N and 81S (provided the width of the map does not 
exceed 28 1/2"), MEGRI also sets up tables for the dimenstoned varlable 
B(J) which Is used to convert any latitude to Inches for plotting within 
this grid, It Is necessary to call MEGRI before calling MEGRW, ALATI, or 
ALONG. MEGRW draws Inside grids, makes annotations, and draws minute 
marks, ALATI converts a given latitude to Inches for plotting within the 
grid specified by MEGRI, ALONG similarly converts a given longitude to 
Inches, 

The subrout!nes DECID, MECHA, and MEPLO do not require MEGRI to be 
called earlfer, DECID and MECHA are used In subroutines such as MERGR, 
which annotates values along a track, DECID ts used to decide which side of 
the track (above, below, left, or right) the annotation should be made for 
maximum clarity In the annotations, MECHA makes the actual annotatlons, 
MEPLO ts sImilar to the IBM subroutine FPLOT; It makes a cross at the 


plotted point, 


Subroutines MERGP, MERGR, and TROJ2, 
MERGP !s used In programs MGRA6, METC6, and MGGG2 to plot the value of 
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gravity, depth, or magnetic Intensity along the track on a Mercator chart, 
The value can be plotted at right angles to the track or at a fixed azimuth, 
In addition, there are options which allow the plotting of the track Itself 
and Indicate by crosses on the plot the times at which navigation polnts 
occur. 

MERGR Is used In programs MGRA4, METC8, and MGGG8 to annotate the 
value of gravity, depth, or magnetic Intensity along the track on a Mercator 
chart. 

TROJ2 !s used In programs NAGS!, NATS2, and MGPR2 to project values of 
gravity, depth or magnetic Intensity along a Itne with given azimuth, If 
this I!ne passes through the origin at VLAT, VLON, and the azimuth Is ANGLE 


(see Flg, 3) 


(VLAT, VLON) 


ALC ULL 
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then the projected distance DIS Is given by 
DIS = DLON x COS(ANGLE) + DLAT «x S/N (ANGLE) 


where = 
DL ASRS Poh cla eee, Lela] 


DLON = (xLON-VLON) x CosC(CxLAT+VvLAT )/2) if DLATAI 


or 

DLON = (KLON ~VLON) x C SIN (XLAT) - S/N (VLAT)) [ CXLAT- vLAT) 
Ts PLAT Se 

zigs In the orlginal track can cause values of DIS to Increase and decrease, 

thus resulting In overlaps In the projected profile, TROJ2 has been pro- 

grammed In such a way as to eliminate any overlapping by accepting only 


those projected distances DIS that are In Increas!tng order, 
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